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ABSTRACT 

The signal-to-noise ratio (SNR) for black hole quasinormal mode sources of 
low-frequency gravitational waves is estimated using a Monte Carlo approach 
that replaces the all-sky average approximation. We consider an eleven dimen- 
sional parameter space that includes both source and detector parameters. We 
find that in the black-hole mass range M ~ 4-7 x IO^Mq the SNR is significantly 
higher than the SNR for the all-sky average case, as a result of the variation of 
the spin parameter of the sources. This increased SNR may translate to a higher 
event rate for the Laser Interferometer Space Antenna (LISA). We also study 
the directional dependence of the SNR, show at which directions in the sky LISA 
will have greater response, and identify the LISA blind spots. 

Subject headings: gravitational waves — black hole physics 

1. Introduction 

An important source of gravitational waves in the frequency range 10~^-10~^IIz, which 
is the frequency band for the planned space-bourne gravitational-wave detector Laser In- 
terferometer Space Antenna (LISA), is the inspiral and merger of two supermassive black 
holes. The emitted gravitational waveform includes three major parts, corresponding to 
the inspiral, merger, and ringdown of the system. The last stage, the ringdown, is the set- 
tling down of the resulting black hole to quiescence as required by the "no-hair" theorem, 
and is characterized by (an incomplete) set of complex frequencies that depend only on the 
macroscopic parameters of the black hole, namely its mass M and spin angular momentum 
J (NoUert 1999). The computation of the complex frequency is done using black hole per- 
turbation theory (Chandrasekhar 1983), and its detection would be a smoking gun for the 
source being a black hole, and an accurate means for determining its parameters (Berti et 
al. 2006). 
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Naively one would think that black hole quasinormal modes are not a strong source 
of gravitational waves for LISA. Indeed, the typical (real part of the) frequency for a 
Schwarzschild black hole (for quadrupole radiation 1 = 2 and azimuthal mode m = 2 at 
the longest damping time n = 0) is /^=2,m=2,ra=o ~ 1-2 x lO'^(M0/M)Hz whereas the e- 
folding time for the exponential decay of the amplitude (associated with the imaginary part 
of the complex frequency) is r£=2,m=2,n=o ~ hb{M/ Mq)iis. While for a stellar black hole this 
frequency is several orders of magnitude too high to be detected by LISA, the frequency 
for supermassive black holes falls inside LISA's good sensitivity frequency band. The qual- 
ity factor Q associated with the oscillator is Q := tt/t, or for a Schwarzschild black hole 
(5 ~ 2, i.e., within two oscillations the amplitude of the waves is suppressed by a factor 
of e^'" ~ 535. (While for rotating black holes the quality factor is higher, it is still low 
unless the black hole's spin approaches extremality. Indeed, the quality factor can be ap- 
proximated by Q22Q = 0.7000 + 1.4187(1 -j)"°-^^^° ± 0.88% (Berti et al. 2006). A different 
parameterization of Q220, which however agrees with this one to within 1% was first given by 
Echeverria (1998).) Therefore, a black hole is expected to be a very poor resonator, and the 
waves difficult to detect. Regardless of these naive expectations, Flanagan & Hughes (1998) 
showed that the signal-to-noise ratio (henceforth SNR) from the ringdown stage could be 
comparable to that of the inspiral stage, so that waves emitted during the ringdown stage 
may be an important source for LISA. 

The signal to be detected by LISA comes from sources that span a large parameter 
space, including both source parameters and detector parameters. To find the SNR from a 
large sample of ringdown sources as a function of the black hole's mass one needs to include 
our ignorance of the particular parameters. Specifically, we do not know the parameters 6, (f) 
that are the spherical polar angles that refer to the direction to the source, and ip which 
is the angle between the waves polarization plane and the plane = that determine the 
detector pattern functions F~^'^{6,(j),il>) for the two gravitational wave polarization states. 
Specifically, 

F+(6I, 0, V') = -(1 + cos^ 0) cos 20 cos 2^' - cos 9 sin 20 sin 2ip 

F''{e,(f),ijj) = ^(1+ cos^^*) COS 20 cos 2V' cos 61 sin 20 sin 2V'. 

We also arc ignorant of the source's parameters j3 that are the spherical polar angles in the 
frame of the source, defining the orientation of the source in its local frame (i.e., determining 
the orientation of the black hole's spin axis). Other unknown parameters are the black 
hole's dimensionless specific spin angular momentum j := cJ/{GM'^), the phases of the two 
polarization states 0^^^„, the total radiation efficiency erd M'^c"'^ J^{dE / df)df (i.e., 
the fraction of Mc^ being emitted in gravitational waves) and e"*", the fraction thereof in 
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the + polarization state, in addition to the distance to the source, which we take to be 
the luminosity distance D^. Altogether, we have an eleven dimensional parameter space (3 
dimensions for the detector parameters, and 8 for the source's) (Thorne 1987). 

The currently available method is to make the ignorance of the parameters manifest by 
averaging over some of the parameters (the detector parameters, i.e., 6, (f), and ip, in addition 
to the the two source parameters l and f3) and fixing other parameters (the source parameters 
j, £rd, Dl). Consider the question of what is the SNR for ringdown sources of a 

particular mass. This question is important for detection purposes or for estimation of the 
population statistics. (This question is different from the question of what is the SNR for a 
ringdown source with particular parameters, which we address in what follows.) Specifically, 
as one black hole source may be in one direction in the sky with a certain value for the 
specific spin (and similarly for the other parameters), and another black hole source may be 
in a different direction with a different value of the specific spin (and the other parameters), 
one may assume that all black hole sources have average (or fixed) values for all parameters 
(notice that Bcrti ct al. (2006) did not average over the spin angular momentum j, but rather 
fixed it). Hereafter, we refer to this approach as the "all-sky average" (because it involves 
averaging over the detector parameters, although it involves also parameter fixing and not 
just averaging). In this approach one would argue that a good estimate for the total SNR for 
ringdown sources in a given mass interval would be found by taking the spin of all such sources 
to be the average value, say j — 0.5 if the distribution function is taken to be uniform. The 
all-sky average approach was used recently by Berti et al. (2006) to find the SNR. The all-sky 
average approach clearly is useful as an approximation method to find the SNR, and as we 
show below docs it to very high accuracy. Notice, that our ignorance of the parameters to be 
used is of two kinds: first, there arc observational unknowns, that will become known after 
precise observations are made. Second, there are theoretical unknowns, that would become 
better known (including their distribution functions) after better theoretical modeling and 
understanding of the sources' physics is obtained. We do not distinguish here between the 
two types of unknowns. 

An alternative approach is to integrate over a finite number of sources, and use Monte 
Carlo integration to refiect the ignorance of the parameters. That is, one uses random values 
of the parameters for each source, spanning the parameter space, and then averages over a 
large number of such sources. The range of the parameters and the distribution functions 
for the random sample then reflect the physics of the source and the gravitational-wave 
generation. In practice, we take here for simplicity all distribution functions to be uniform, 
but this assumption is easy to relax if evidence is presented for non-uniform distributions. 
Specifically, we take the sources to be homogeneous and isotropic (so that the uniformity 
of the distribution function is in D\ and cos^), and the sources' spin axes to be isotropic 
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(i.e., the distribution function is homogeneous in cost). Recently, our understanding of 
gravitational-wave sources has improved dramatically. We now know that the radiation 
efficiency for binary inspiral, at least for the variables simulated using numerical relativity, 
is ~ 3% with circular polarization. However, it was argued that Cj-d ^ 7% for some cases 
(Berti et al. 2007), and for head-on collision it is e^-d ~ 0.1% with + polarization (Berti et 
al. 2006). More recent numerical relativity simulations for the merger of comparable mass 
black holes have found a relationship of the gravitational-wave radiation efficiency to the 
spin of the initial black holes, and a relationship of the spin of the final black hole to the spin 
of the initial black holes (CampaneUi et al. 2007). For the specific case of initially ahgned 
or anti-aligned black hole spins, the relationship of the radiation efficiency in gravitational 
waves to the final black hole's spin j is given by 



as obtained from a quadratic fit of the numerical relativity results of CampaneUi et al. (2007). 
The final black hole's spin was found to be in the range 0.35 ^ j ^ 0.95, which implies 
radiation efficiency in the range 2.2% ^ erd ^ 8.3%. The dependence erd = ^rd{j) reduces 
our parameter space from 11 to 10 independent dimensions. In what follows, we specialize 
our discussion to the case of quasi-normal modes from black hole sources resulting from the 
merger of black hole binaries with initially aligned or anti-aligned spins. This specialization 
allows us to use actual numerical relativity values for the radiation efficiency in gravitational 
waves. Our method can be applied also for black holes resulting from less restrictive initial 
spin configurations when corresponding numerical relativity results become available. 

One may naively expect that in the limit of very many sources, the Monte Carlo results 
reduce to the all-sky average results. However, they do not. The large number limit of the 
Monte Carlo integration approaches the all-sky average only if the parameters of the SNR 
are independent. When considering only Schwarzschild black holes, or even all black holes 
with a fixed value for the specific spin angular momentum, indeed all the parameters are 
independent, i.e., they live in different spaces. However, when one considers that for any 
given value for the black hole mass M the specific spin spans the range ^ j ^ 1 (or, more 
specifically for the case considered here in detail, 0.35 ^ j ^ 0.95), these parameters are no 
longer independent. 

Specifically, the SNR squared is given by (Berti et al. 2006) 



(j) = 0.994 - 0.438 ^4-751 - 3.74 j - 0.4545 j , 
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Here, b± = (l/Tfmn)/[(l/''"^mn)^ + 47r^(/±/£m„)^] is the Fourier transform of the time domain 
single frequency quasi-normal signal that following Berti et al. (2006) we approximate using 
the delta-function (in frequency) approximation, and Sh{f) is the LISA noise spectral den- 
sity. In practice we calculate Sh{f) using the approximation given by Finn & Thorne (2000) 
and by Larson et al. (2000), which agrees well, except for oscillations at the high-frequency 
end, with the LISA Sensitivity Curve Generator resultil]. The amplitudes in the two polariza- 
tion states are denoted by Af^^, and Sirnnii^, (3) are the spin-weighted spheroidal harmonics 
of spin weight 2, which depend on j through the overtone index n. Also the amplitudes Af^^ 
are functions of j, as 



and all three variables /^^n, Qimn, and e^d are functions of j. Clearly, (^^^^(j) Semn{j)) ^ 
(Aimnij)) X {SimnU)) whcu the averaging is done over j. We may not therefore expect the 
qualitative dependence of as a function of the mass M to remain unchanged when the 
all-sky averaging assumption is relaxed. In this Paper we study this question, and find the 
SNR when the sources are obtained from a random sample. That is, we do not assume 
following Berti et al. (2006), that all black hole quasi- normal mode sources have the average 
or fixed values for all parameters, and we let these parameters take random values, as one 
would expect from an actual sample of sources for LISA. The all-sky average approximation 
of Berti et al. (2006) assumed that = = i, F+F^ = {F+F"") = 0, and 

that S^^^ = Sg^^Simn = {Simn ^imn) = 4^, wherc an asterisk denotes complex conjugation. 
It was further assumed following Flanagan & Hughes (1998) that = (p^rnn ~ 

that the amplitudes = A^^^ = A^rnn- The specific spin angular momentum j and 

the radiation efficiency e^^ were fixed. To compare our Monte Carlo results results with the 
all-sky average approach, we set the spin angular momentum in the all-sky average case to 
the value that gives the average radiation efficiency. Specifically, we find j that would solve 
erd(j) = (erd(j)) = Sj2in ^i-dO') Omax - jmin)- Ovcr the range of j for which erd(j) is valid, 
we find in practice that (erd(j)) = 0.03824 and the corresponding spin angular momentum 
is j = 0.71882. 

Our Monte Carlo approach is a simple application of Monte Carlo integration, as we 
are using uniform distribution functions for our random variables (see below for justifica- 
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-'^The LISA Sensitivity Curve Generator (SCG) is available online at this URL: 
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tion). The (pseudo-)random number generator is based on the Mersenne Twister algorithm 
(Matsumoto & Nishimura 1998). While not considered secure for cryptography applications, 
this generator is certainly reliable enough for our purposes. Specifically, we sample points 
uniformly from the integration region to estimate the integral and its error. Namely, we 
calculate the SNR by estimating the integral in ([2]) by summing the integrand over the ran- 
domly chosen sample, dividing by the volume of the parameter space and by the number 
of iterations. The estimated error cr in the Monte Carlo integration is done by calculating 
the variance, according to cr^ = (V/N)'^ Yl!i=i[p{^i) ~ (p)]^5 where V is the volume of the 
parameter space and N is the number of elements in the sample. Here Xi denotes collectively 
the eleven parameters for the SNR for the i^^ member of the sample. The reason why this 
very simple approach to Monte Carlo integration is sufficient for our needs is that we choose 
the sources to be uniformly distributed across the parameter space. Specifically, we take the 
sources to be distributed homogeneously and isotropically, which seems to be a very reason- 
able assumption due to the cosmological nature of the sources. We also take the spin axis 
of the black hole members of the sample to be distributed isotropically. These assumptions 
mean that we take D\, cos 6', cost, and (3 to be uniformly distributed. We have little 
motivation to prefer a non-uniform distribution function for the other parameters. Specifi- 
cally, we currently do not have enough numerical relativity results to determine a realistic 
distribution function for the parameters ?/',<^"'"'^ or in addition to the population statis- 
tics of j. We therefore take these parameters to be distributed uniformly in their respective 
parameter spaces. The radiation efficiency erdO) is determined according to Eq. ([1]). 

Our approach allows us also to study sources with particular parameters instead of the 
all-sky averaged counterparts. Specifically, we consider the directional variations of fixed 
sources, and show in which directions in the sky LISA would be more sensitive, including 
finding LISA blind spots. 

In all our simulations we take a standard zero-curvature cosmological model with the 
5-year WMAP values VLm = 0.279, = 0.721, and Ho = 70.1 kms~^Mpc"^ (Hinshaw et 
al. 2009), and the luminosity distance - redshift relation to be given by (Carroll et al. 1992) 



2. Testing the Monte Carlo integration in the Schwarzschild and Kerr cases 

We first test our Monte Carlo integration for the simple case in which all black holes are 
assumed to be Schwarzschild black holes (j = 0). In particular, we confront our results with 
the all-sky average results, and study how well they agree. Figure [T]A shows the SNR for the 
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all-sky average and for the Monte Carlo simulation. We find that the two signals overlap, 




Fig. 1. — Upper panel (A): The SNR for the all-sky average (results analogous to those 
of Berti et al. (2006)) (dashed curve) and the Monte Carlo results with = 10 000 
Schwarzschild sources (solid curves and error bars) at luminosity distance of 3Gpc (red- 
shift of 2; = 0.522). Radiation efficiency for the all-sky average case is 3%. Lower panel (B): 
The SNR for the all-sky average (results analogous to those of Berti et al. (2006)) (dashed 
curve) and the Monte Carlo results with = 1 600 Kerr sources (solid curve and error bars) 
at luminosity distance of 3Gpc (redshift of 2; = 0.522), for j = 0.8. Radiation efficiency for 
the all-sky average case is 4.948%. 

as indeed is expected. In particular, in the Schwarzschild case our Monte Carlo results agree 
to the expected accuracy level with the results of Berti et al. (2006). 

The Schwarzschild test does not test all sectors of our code, as it does not calculate the 
spin-weighted spheroidal harmonics. Fixing the value of j we can test the remaining sectors 
of the code. Figure [TJ3 displays the SNR for the all-sky average (for a fixed value of j) and 
the Monte Carlo results, which are in agreement to within the error bars. 

We also test the convergence of our Monte Carlo code, letting all eleven variables be 
random. Figure [2] shows the SNR for different Monte Carlo runs, with different sample sizes. 
The convergence of the result is indeed as expected, namely scaling with A/iV, A^ being the 
number of sources in the sample (see a more detailed analysis of Fig. [2] in Section [3]). 
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Fig. 2. — The SNR for the Monte Carlo results with Kerr sources at variable luminosity 
distances out to redshift z = 2 and variable y. N = 800 (thin solid curve) and N = 1600 
(thick solid curve). The corresponding all-sky average results are shown for comparison, at 
redshift z = 1.509 (dashed curve). The insert shows the ratio of the Monte Carlo to the 
all-sky average SNR. 



As we already pointed out, we do not expect the behavior of the SNR when j is let 
to vary to necessarily be identical to the case of j being averaged. In Figure [3] we fix all 
the variables, and let only j be a random variable. The distance of the sources is fixed at 
redshift of z = 1.509. Indeed, we find that the behavior is not identical. While for very low 
and very high masses the behavior is nearly indistinguishable, there is a range of masses, 
specifically M ~ 4-7 x lO^M© for which the SNR in the variable case is significantly higher 
than in the all-sky average case. The ratio of the two SNRs spikes at M ~ 5.6 x lO'^M©, 
where p/paii-sky ~ 1-4 which would translate to an event rate greater by a factor of ~ 2.75. 
This increase in the SNR is 14.2ct, so that it is not an artifact of the Monte Carlo error 
associated with the finite sample size. 

Next, we let all the variables — except for the distance, which we fix at redshift of 
z = 1.509 — be random. Figure H] shows the SNR, which preserves the behavior we found in 
Fig. [3l Specifically, the sharp increase in SNR for masses in the range M ~ 4-7 x lO^M© 
is retained. The SNR is increased by a factor of 1.6, that translates to an increased event 
rate by a factor of 4.1. This increase is 7.4(j. Indeed, this is the expected behavior, as all 
other variables live in independent spaces. Notice, that as the number of dimensions here is 
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Fig. 3. — The SNR for the all-sky average (dashed curve) and the Monte Carlo results with 
variable j for Kerr sources for fixed parameters other than j: N = 800 (thin solid curve) 
and = 1600 (thick solid curve). The insert shows the ratio of the SNRs. The all-sky 
average results are obtained for j = 0.71882. Here, all sources are at redshift z = 1.509, 
COS6 = cose = 0.7795, P = (p = 0.1894, = 0, 0+ = 0^ = 0, e+ = e^ and e^d = 3.824%. 

higher than in Fig. [3l the errors are commensurately larger. 

We present in Fig. [2] the results for all variable being random, including the distance, 
from the distance to the Virgo cluster (18.0 ± 1.2 Mpc) up to redshift of 2; = 2. We find 
that the SNR is greatly increased compared with all sources being at the average distance 
(squared) (redshift of z = 1.509) over a large range of masses, specifically for M ^ 3.5 x 
IO^Mq. The increased SNR peaks at M ~ 5 x IO^Mq, with an increase in the SNR by 
a factor of 2.4, which translates to an increase in the event rate by a factor of 13.8. This 
increase is 9. la. As LISA will detect sources at all distances, we believe the increased SNR 
shown in Fig. [2] is closer to the actual SNR for LISA. 

Most importantly, our approach allows us to find the SNR in different directions in 
the sky. In Fig. [S] we show the SNR as a function of the colatitude 6 for various values 
of (j), along with the all-sky counterpart. We show that the same source (namely, same 
source's parameters) located in different direction have very different SNR projected on the 
LISA detector. In particular, there is a LISA blind spot at 6 = 7t/2 and (p = it/A (and 
other values of on the LISA equatorial plane given by the rotational symmetry, namely at 
= 37r/4, 57r/4, Tn/A; our numerical results are within 0.1% of these values.). 

This work can be extended in the following ways. First, one may relax the assumption 
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Fig. 4. — The SNR for the all-sky average (dashed curve) and the Monte Carlo results with 
all variables random (except for distance that is fixed at redshift of z = 1.509) for 800 Kerr 
sources (thin solid curve and error bars) and for 1 600 Kerr sources (thick solid curve and 
error bars). The all-sky average results are obtained for j = 0.71882. 

of uniform distribution functions for the random variables, and base the variables on more 
reahstic distribution functions. Second, we use here the delta-function approximation to 
the quasinormal mode frequency. This assumption may be relaxed (Berti et al. 2006), but 
the change in the SNR is not expected to be significant, and third, one may extend the 
single-mode analysis done here to multiple modes (Berti et al. 2006). 
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